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Abstract: The simple entanglement of TV-body TV-particle pure states is extended to the more general 
M-body or M-body TV-particle states where TV ^ M. Some new features of the M-body TV-particle pure 
states are discussed. An application of the measure to quantify quantum correlations in a Bose-Einstien 
condensate model is demonstrated. 
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^ ■ 1 Introduction 

Entanglement plays an important role in quantum information theory and quantum computation [l]-[3]. 

O ■ 

\Q , Entanglement normally considered is about TV-body TV-particle states. The most famous and widely 
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studied of these are probably the Bell states [4] with TV = 2, and TV-body TV-particle GHZ [5] and W [6] 
states with TV > 3. Often, however, that the number of particles TV is not equal to the number of kinds 
of particles M in a quantum state. Such multipartite entangled states can be generated naturally in 
numerous many-body quantum systems, such as spin systems [7]- [9], one dimensional chains or lattices of 
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various dimensions of Bose or Fermi many-body systems [10]-[12] including two-body TV-particle entangled 
states generated in the Bose-Einstein condensate [13] with sufficiently large TV, of which physical systems 
can be fabricated with modern technology. Specifically, entangled states can be realized by burning the 



system into a desired configuration and then cooling it down, which leads to the so-called ground state 
entanglement [14]. In order to classify different entangled states and study the entanglement of these 
types of states, it is necessary to extend the definition of entanglement to M-body TV-particle pure states 
where TV ^ M. Up till now, there have been at least four different approaches to the problem, namely, 
the entanglement of modes [15], the quantum correlations [16], the entanglement of particles [17], and 
entanglement witness [18]. The entanglement of modes proposed by Zanardi [15] is defined in terms 
of local entanglement in terms of local reduced density matrix in Fock representation. The quantum 
correlations for pure states introduced in [16] are determined by the so-called Slater rank of the state. 
The entanglement of particles provides an operational definition of entanglement between two parties who 
share an arbitrary pure state with a standard measure [17]. Finally, entanglement of witness is defined 
in [18] based on witness operators. Each of these approaches emphasizes one or several respects of the 
problem. There is still no consensus reached on which one is the most suitable. 
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As an extension of the entanglement measure for iV-body iV-particle pure states [19], in this paper, 
we propose an entanglement measure for M-body iV-particle states for any M and N, which can be 
regarded as an extension of the entanglement of modes proposed by Zanardi [15], but does not count 
in entanglement among local identical (indistinguishable) particles. In Sec. 2, we will propose the 
entanglement measure for M-body iV-particle pure states for any M and N. Some simple examples 
taken from quantum many-body systems will be analyzed with this measure in Sec. 3. Sec. 4 is a short 
discussion. 

2 The entanglement measure for M-body TV-particle pure states 

For a system consisting of M kinds of particles, any pure state \*&) can be expanded in terms of basis 
vectors \n\, o\\ nz, 02; • • • ; Hm, &m) in the tensor product space Vk t <8> Vk 2 ® • • • ® Vk M as 

l*>= E C?™:::™\n 1 ,v 1] n 2 ,*2 ] --- ] n M ,a M ), (1) 

Ill(7lIl2(T2-llMITM 

where rn (1 < i < M) is the number of the z-th kind of particle, Oi = 0, 1, • • • , fcj for 1 < i < M denotes 
the internal degrees of freedom, and C™ 1 "*."'"^ is the normalized expansion coefheient. Then (1) will be 
called an M-body pure state. In addition, if the total number of particles N — Y^iLi ™i ^ s a conserved 
quantity with 

N\V)=N\V), (2) 

(1) will be called an M-body A^-particle pure state, which will be denoted as |\P(M, N)) in the following. 
According to the above definition, Bell states are two-body two-particle states, GHZ and W states are 
iV-body A^-particle states with N > 3. The corresponding density matrix of (1) is 

P* = |*><*|. (3) 
Let a\ a (a ia ) be local particle creation (annihilation) operators that satisfy 

[a.ia,a) a ,]± = a ia a} ja , ± a] a ,a la = 5ij5 aa >, (4) 

[out, a ja >]± = [at,, a] a ,]± = 0, (5) 
for fermions or bosons. The wavefunction \$>{M,N)) can be expressed as 

\9(M,N)) = e c niai ... nMa M)?M)%---(«M)*z\Q), (6) 

where 10) is the vacuum state and (aj)™ 4 = Y]„i„i „ t Q a \ < , , at t a}. 4 • • • 4 , in which B a \ t f is 
the expansion coefficient. 
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In (1) and (6), {\rii<Ji}} is recognized as a set of the local basis vectors for the i-th kind of particles, 
which can be used, for example, to describe a set of local states on the i-th site of a lattice in many 
lattice models. In such cases, only entanglement among particles on different sites is of interest, which 
can also be used to signify quantum correlations in a model system. Therefore, entanglement among 
particles on the same site will not be considered. Similar problems were also considered in many works. 
For example, hierarchic classification for arbitrary multi-qubit mixed states based on the separability 
properties of certain partitions was considered by Diir and Cirac in [20]. 

Under the replacement at 4 — > X; « , where X { » is simply a symbol, the operator form in front of the 

ta j 3 3 

vacuum state on the right-hand-side of (6) becomes a homogeneous polynomial of degree N in terms of 
the {XJ, 

F C (X 1 ,...,X M )= J2 C ai ... a „{X 1 )%---(X M )™. (7) 

(71---CTJV 

It should be understood that Xj is a multi- value symbol with Xj = X irT i for any a). An alternative 

3 J 

definition of entangled states with respect to the local bases can be stated as follows: The state \^(M, A)) 
is an M-body A-partite entangled state if the corresponding polynomial Fc{X-i, • • • , Xjf) on complex 
field C cannot be factorized into the following form 

Fc(X-i, ■ ■ ■ , X M ) = Fa(^% 1 , • • • , Xj ! )i ?1 B(Xj i+1 , • • • , Xj M ) (8) 

for 1 < I < M — 1, where {ii ^ i 2 ^ ■ ■ ■ ^ im} can be in any ordering of {1, 2, • • • , M}. Otherwise the 
state |^(M, N)) is not an M-body A/'-particle entangled state with respect to the M sets of local bases. 
The state |^(M, N)) given in (6) is disentangled (separable) if the polynomial Fc can be factorized into 
a product of polynomials of Xj as n^ii ^Ai(Xj). In other cases, the state is partially entangled. 

As is shown in [19], a criterion for distinguishing whether a homogeneous polynomial is factorizable 
can be established by using the von Neumann entropy of the reduced local density matrices. As an 
extension of the entanglement measure [19] for A^-body A-partical qubit system, we define entanglement 
measure for the M-body pure state (1) as 

f wEf if Si ?0 V i, 

rjM = (9) 
t otherwise, 
where M is the number of different kinds of particles involved in (1), 

S {i) =-Tr[p (i) log Pi (p (i) )] (10) 

is the extended local reduced von Neumann entropy expressed in terms of local reduced density matrix 
for the i-th kind of particles only obtained by taking the partial trace over the subsystem, and pi is 
the total number of the Fock states of the i-th kind of particles in the system. We use the logarithm to 
the base pi instead of base 2 used in qubit system [19] to ensure that the maximal measure is normalized 
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to 1. Hence, the maximal value of the measure is 1 regardless how many particles and how many kinds 
of particles are involved in the system. While the measure is zero when a pure state is not a genuine 
M-body entangled state. 

Similar to the JV-body TV-particle case [19], the measure (9) quantifies genuine M-body entanglement. 
The measure is zero if Si = for any 1 < i < M, in which the corresponding state is not a genuine 
M-body entangled one. In contrast to the measure defined in terms of local von Neumann entropy [16], 
which provides information of local entanglement only, the measure (9) provides information about overall 
quantum correlations among M kinds of particles in a quantum many-body system. 

In the pure state (6), the local particle number is not a conserved quantity in general. In such a case, 
the local reduced density matrix pu\ is built in the local Fock space spanned by {\rnai}} (0 < rii < pi) 
with the following block-diagonal form: 



P(i) = 



/(ni = 0) \ 

(m = i) 

(m = 2) 



(ii) 



V {n l =p l )J 
where each sub-matrix (rij = q) is spanned by the g-particle Fock states with the internal degrees of 
freedoms < cti < fcj. Therefore, the extended reduced von Neumann entropy Su\ is not only invariant 
under the Local Unitary (LU) transformations for each block in (11) with respect to internal degrees of 
freedoms {uj} for fixed rij, but also invariant under unitary transformations for the subspace spanned 
by the entire Fock states with the local particle numbers < rii < pi. Hence, the measure (9) is also 
invariant under both the Local Unitary (LU) transformations with respect to internal degrees of freedoms 
and unitary transformations for the subspace spanned by the entire Fock states. In addition, the measure 
does not count in entanglement among local identical particles with the same label i and different internal 
degrees of freedom because the local basis vectors are recognized to be {\rii(Ji)}. 

In the M = N qubit cases discussed in [19], the LU invariance is equivalent to the invariance under 
Local Operations assisted by Classical Communications (LOCCs). In such case the measure (9) is both 
LU and LOCC invariant [21]. In N ^ M cases, however, LOCCs may be useless in any quantum 
information processing protocol. For example, LOCCs will not only keep the entanglement the same, 
but also be irrelevant to any quantum information processing processes for M ^ N. In addition, though 
(9) is invariant under local unitary transformations with respect to different number of particles, such 
transformations will result in local state mixing among different number of particles, which obviously 
violates the super-selection rule, and thus is impossible [22, 23, 24]. However, the measure (9) is indeed 
useful in quantifying the M-body TV-particle entanglement, of which some examples will be shown in the 
next section. 
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3 Some simple examples 



Using the above definition, one can easily write out many non-trivial states with M ^ N which are 
quite different from those with M — N . For simplicity, let us assume that there are M kinds of spinlcss 
fermions or bosons. Then, the simplest case is a M-body one-particle state with 

MM M 

|*(M,1)) ^C 1 \1Q^}+C 2 \oTo^) + --- + Cm\0^^), (12) 

where |0 ■ ■ ■ 010 • • • 0) with 1 in the i-th position denotes one-particle state for the i-th kind of particle 
and vacuum state of the remaining particles, and Cj (i = 1, 2, • • • , M) are normalized non-zero expansion 
coefficients. These types of states emerge naturally from many models, such the Bose- and Fermi-Hubbard 
models. 

In TV-body A-particle cases, generally a state is entangled with respect to its intrinsic degrees of 
freedom, such as the third component of spin of photons or atoms. In such cases, the measure (9) is 
invariant under LOCCs or LUs. In N ^ M cases, however, LOCCs may be useless in any quantum 
information processing protocol. In order to make this point clear, for example, let us replace spinlcss 
fermions or bosons in (12) by those with internal degrees of freedom <jj = 1, 2, • • • , fa for i = 1, 2, • • • , M. 
Then, (12) becomes 

MM M 

\^>(M, l)(<7i<72 • • • <7 M )) -Ci|1^0^+C 2 |0l7l0~^ + --- + CAf|^0O- (13) 

fa should be all the same for any i if these are identical particles, otherwise they are not identical. In this 
case, any LOCC between i-th and j-th parties is equivalent to a corresponding LU transformation with 
respect to the internal degrees of freedom labeled by <ii or aj . After any such LU transformation, say 
at i-th party, the vacuum state of the i-th particle for any i remains invariant, while local one-particle 
state still remains to be one-particle state with a linear combination of different internal labels. The j-th 
party can not get any information after local operation at i-th party assisted by classical communication 
between them and then measuring the resultant j-th local state because the j-th state remains to be the 
same after any local operation at i-th party. Therefore, LOCCs will not only keep the entanglement the 
same for the state shown in (12), but also be irrelevant to any quantum information processing processes 
used in M = N cases even if those M Fock states have additional intrinsic degrees of freedom as shown 
in (13). Actually, N ^ M entangled states, such as that shown in (13), are entangled with respect to 
local particle numbers in M different local Fock states and not with respect to local intrinsic degrees 
of freedom. For a state, such as that shown in (13), effective local operations similar to those used 
in M — N cases are local particle number non-conserved unitary operations, namely local projective 
operations, a\0) + b\l), where |0) stands for vacuum state, |1) stands for one particle state, and a, b € C 
with |a| 2 + \b\ 2 — 1. As for entanglement in spin degrees of freedom, in which any local operation may 
violate conservation of spin in the third component, any local particle number non-conserved unitary 
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operation violates both local and total particle number conservation. Such operations obviously violate 
the super-selection rule, and thus are impossible to be implemented. While the entanglement measure 
(9) is invariant under any such local particle number non-conserved operation. 

Although local particle number non-conserved operations are physically impossible, measure (9) can 
be used to quantify entanglement of any M-body entangled state. For example, according to (9), the 
entanglement measure for (12) is 

M 

V (M > N=1) = ~ £ W lo S 2 (l^| 2 ) + (1 - l^| 2 )log 2 (l - |a| 2 )) . (14) 

i=l 

It can be verified by conditional maximization that (12) reaches its maximal value when |C$| = 1/yM 
with M > 2. Clearly, in the cases with maximal measure 

^ JV=1) =log 2 (M)-^^log 2 (M-l) for M>2. (15) 

(2 1) 

It is similar to the Bell states with ?7max = 1 when M = 2, and is equivalent to the W states with 
?7max = 0.918296 when M = 3, which is the same value as that of a W state for a qubit system calculated 
in [25]. 

Though there is controversy about whether such states as that shown in (12) are entangled or not [26, 
27], it was clearly demonstrated in [28] that such state is indeed a single-particle entangled state, which 
was also used in a two-party protocol for quantum gambling [29], and is proved to be useful in the 
so-called data hiding protocols [24] . The measure (9) is in agreement with those observations. 

Another example of our application of the measure (9) is to study quantum phase transition in a 
Bose-Einstien condensate model with n bosons in a two-dimensional harmonic trap interacting via a 
fairly universal contact interaction [30, 31, 32], of which the Hamiltonian can be written as 

oo 

H = y^wja]aj + g £ ^+j\fc+iV r i j;fc,za|aja fe a;, (16) 

j=0 i,j,k,l 

where a\ (at) is creation (annihilation) operator of the i-th boson, and Vij-ki = (k + l)\ / (i\j\k\VA k+l )^ . 
In this model, the total number of bosons h and the angular momentum L with 

oo oo 

n = ^a\ ai , L = J2ja]a 3 (17) 

i=a j=a 

are two conserved quantities. Therefore, for given n and L, eigenstates of (16) can be written as 

\n,L;r)= £ ^^^^^ 

n 2 n 3 ---n L 

a\ L -^^4 n *at n ° ...a^\0), (18) 

where c^ 2 „ ... are expansion coefficients, and r is an additional quantum number needed in distinguish 
different eigenstates with the same n and L. 
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To explore transitional patterns in the system, we calculate the entanglement according to (9) along 
the yrast line up to L = 6, in which the local reduced density matrices pu\ for the i-th kind of bosons 
are obtained by taking the partial trace over the subsystem. Fig. 1 shows the entanglement measure (9) 
of the system as a function of the total number of bosons n along the yrast line with L = 2, 3, • • • , 6. 
Generally speaking, the maximal value of the measure decreases with L; and the measure for any L 
decreases with increasing of n due to the s-boson dominance in the yrast states when n is sufficiently 
large. It is obvious that there is always a peak in the measure at or near n = L + 1 when L > 4. The peak 
in the entanglement measure signify that there is a quantum phase transition occurring near n = L + 1, 
in which the total number of bosons n serves as the control parameter. 

In order to show that the peak point in the entanglement measure shown in Fig. 1 indeed is a critical 
point of the system, we study yrast state occupation probabilities £j (n) for bosons with different angular 
momentum I as a function of n according to 

<j>i{n) = {a\ai)/n. (19) 

The results for L = 2, 3, 4, 5 cases are shown in Fig. 2, which shows that the values of the occupation 
probabilities for bosons with different angular momentum vary drastically with increasing of the total 
number of bosons, especially for s- and p-bosons. When n = L, the /-boson, especially the p-boson 
components contribute significantly to the yrast states. There are cross points for yrast state occupation 
probabilities 4>i(n) with different I, which are indeed all near n = L + 1 point. The L = 4 and L = 5 
cases shown in Fig. 2 are consistent to the critical points at the entanglement measures shown in Fig. 1 
though it is not shown by the entanglement measures in Fig. 1 for L = 2 and 3 cases. The results show 
that there arc two different phases. One is the p-boson dominant phase, while another is the s-boson 
dominant phase, which are controlled by the total number of bosons in the system. Therefore, the peak 
in the entanglement measure indeed signifies that there is a quantum phase transition occurring near 
n = L + 1, in which the total number of bosons n serves as the control parameter. The result also 
shows that the measure (9) is indeed suitable to benchmark the quantum correlations among bosons with 
different angular momenta. 



4 Discussion 

In conclusion, entanglement of M-body and M-body A^-particle pure states is defined. The simple 
entanglement measure formerly defined for M — N multipartite pure states is extended to the general 
cases with any M and N. The feature of the measure is that it only count in entanglement among 
non-identical particles. For example, the measure does not count in entanglement among particles with 
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Figure 1: The entanglement measure as functions of the total number of bosons n along the yrast 

line up to L = 6. 




Figure 2: Occupation probabilities 4>i(n) with I = 0, 1, • • • , L as functions of n along the yrast line up to 
L = 5, where the upper left panel is for L = 2, the upper right one is for L = 3, the lower left one is for 
L = 4, and the lower right one is for L — 5, respectively. 



different internal degrees of freedom on the same site in lattice models, and so on, which seems to be a 
natural way to define entanglement in quantum many-body problems. 

Though the measure is defined for pure states, it can easily be extended to the measure for mixed 
states p by the convex roof construction with [33] 

rj(p)= mm J2qiV(\*i)(*i\), (20) 

i 

where 3> is the set of ensembles realizing the density matrix p with $ = ({^, \^i}} \ p — Qi\^%}{^i\, 

An example of application of the M-body n-particle pure state entanglement measure to the Bose- 
Einstien condensate model with n bosons in a two-dimensional harmonic trap interacting via a fairly 
universal contact interaction is provided, which shows that the measure is effective to quantify entangle- 
ment or quantum correlations among different kinds of particles in the system. Therefore, the proposed 
entanglement measure is useful not only in quantifying entanglement, but also in studying quantum phase 
transitions in quantum many-body systems with sufficiently large number of particles. 
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